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Abstract 



We investigate the feasibility of the method for reconstructing the equation 
of state and the effective potential of the quintessence field from SNIa data. 
We introduce a useful functional form to fit the luminosity distance with good 
accuracy (the relative error is less than 0.1%). We assess the ambiguity in re- 
constructing the equation of state and the effective potential which originates 
from the uncertainty in Qm- We find that the equation of state is sensitive 
to the assumed VLm-, while the shape of the effective potential is not. We also 
demonstrate the actual reconstruction procedure using the data created by 
Monte-Carlo simulation. Future high precision measurements of distances to 
thousands of SNIa could reveal the shape of the quintessential potential. 

PACS numbers: 98.80.Es; 98.80.Cq 
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I. INTRODUCTION 



Recent various observations 0,@], in particular distance measurements to SNIa HU, 
strongly suggest that the universe is currently dominated by a positive vacuum energy 
density with negative pressure. The smallness of the vacuum energy density ~ (10~ 12 GeV) 
has revived the idea that the cosmological "constant" is not really a constant but rather 
decaying. The idea of quintessence || (see also and references therein) is that vacuum 
energy density is played by a scalar field rolling down the almost flat potential similar to 
cosmological inflation, and a lot of models have been proposed so far. However, there is 
currently no clear guidance from particle physics as to which quintessence models may be 
suitable. Then it should be the observations that decide which model is correct or not. As 
a bottom-up approach, we have proposed that distance measurements to SNIa may allow 
one to reconstruct the equation of state of the dark energy or the effective potential of the 
quintessence field [|J. 

Future observational project, such as SNAP(SuperNova/Acceleration Probe)[| could 
gather 2,000 SNIa in a single year and could put significant constraint on the cosmolog- 
ical parameters (including the equation of state of the dark energy w). In view of the future 
prospect for high-z SNIa search, we investigate in detail the feasibility of the method for 
reconstructing the equation of state and the effective potential of the quintessence field from 
SNIa data 



II. PARAMETERIZING THE LUMINOSITY DISTANCE 

Considering the future prospect for high-z SNIa search, we believe it particularly useful 
to fit the observed luminosity distance di(z) to a function of z which has the following 
properties: (l)the good convergence (the relative error is hopefully less than 0.1% because 
the distance error expected from SNAP will be less than a few percent) for < z < 10; (2) 
the correct asymptotic behavior for z ^> 1 (H(z) oc (1 + z) 3 / 2 ). We present a fitting function 
for d L (z). 

We restrict ourselves to a flat FRW universe henceforth and assume Einstein gravity^ 
In a flat model, the luminosity distance di(z) is written in terms of the coordinate distance 
r(z) to an object at z as 

d L (z) _ , fa dtf r z dz' r 1 2dy' 
JT+zj ~ T{Z} ~ It a!f)~ Jo IMY) ~ J y H(y')y'^ ( 

where to is the present time and y = 1/^/1 + z. It is interesting to note that in terms of y, 
r(y) is a linear function of y if Qm — 1- Therefore, the d 2 r(y)/dy 2 contains the information 



1 http:/ /snap. lbl.gov 

2 While our work was being completed, we became aware of related work [11| which focuses on the 
uncertainty in the reconstruction of the equation of state of dark energy. 



3 The reconstruction equation for the so called extended quintessence [ 12 1 is presented in [13]. 
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of the non-zero pressure (see Eq.@ below). We will elaborate on the advantage of using y 
and r in a separate paper. 

In analogy with Pen's powerful fitting formula for r(z) for a flat FRW universe with a 
cosmological constant H14H , we propose to fit r(z) in the following functional formf]: 



H r(z) = 77(1) -T](y) 
V(y) = 2a 



y- 8 + Py~ s + 7?r 4 + by' 2 + a 



(2) 
(3) 



The requirement H(z)/Hq — > 1 for z — > imply that a is found to be a dependent parameter: 
a = (a (1 + 3/5/4 + 7/2 + 5/4)) 8//9 — 1- /3 — ^ — 5. However, we shall treat a as a free 
parameter for simplicity in numerical calculations. For y — > (z — > 00), we have H(z)/H — > 
l/(a|/ 3 ). For the Einstein-de Sitter universe, a = l,/3 = 7 = 5 = 0. 



A. demonstrating the goodness of fit 

We shall demonstrate that the fitting function Eq.([|) indeed does a good job. For this 
purpose, we calculate the maximum relative error in r(y) between the actual value and that 
calculated from the fitting function Eq.(|3|) in the range 0.3 < y < 1 (10 > z > 0). We 
consider cosmological models consisting of matter and dark energy of constant equation of 
state w = px/px = 0,-1/3,-2/3,-1 with 0.1 < Qm < 1- We fit each template luminosity 
distance, r,, by the functional form Eq.([3|) using the Davidon-Fletcher-Powell method. We 
minimize X^( r i — r {Ui)) 2 ' l r \ with N=30 data points. The result is shown in Fig. 1. We 
find the maximum relative error is less than 0.05%. Our fitting function thus seems more 



powerful by order of magnitude than the one proposed in [10 . 



4 An extension to an open or closed model is immediate once Qk is known: Hqv(z) = 
^fc| _1//2 sinx (jf^l 1 / 2 ^!) — rj(y)fj, where sin^(x) = sin(x)(sinh(x)) if K = 1(— 1). 
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FIGURES 




FIG. 1. The maximum deviation in the range 0.3 < y < 1 (10 > z > 0) between the actual 
value and that calculated from our fitting function is shown as a function of Qm- 



III. RECONSTRUCTING THE EQUATION OF STATE AND THE EFFECTIVE 

POTENTIAL 

One of the prime interest in the reconstruction issue is whether the effective equation of 
state of the x-component w = Px/px is different from —1, which is the unique signature 
of the dynamical vacuum energy. If observations would suggest w ^ — 1, then the next 
urgent project would be the real reconstruction of the effective potential, which should have 
profound implications for particle physics as well as for cosmology. 

In this section, we reconstruct the equation of state of dark energy and the effective 
potential of the quintessence field. Regarding the energy density of the Universe, we as- 
sume two components: nonrelativistic matter with its present density parameter Qm and 
x-component with fix = 1 — ^m- An extension to include the radiation component is 
immediate but with negligible effect. 



A. reconstructing the equation of state 

In terms of a dimensionless variable f = H r, the equation-of-state of x-component w is 
written as 

3dr(z)/dz + 2(1 + z)d 2 r(z)/dz 2 

w(z) = ^7 o \ ( 4 ) 

3 {dr{z)/dz) (Q M (dr(z) dzf (1 + zf - l) 
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—Ayd 2 ?/ dy 2 



3 (df/dy) (n M (df/dy) 2 - 4 



(5) 



w thus depends on the second derivative of the luminosity function. Put another way, the 
luminosity distance depends on w through a multiple integral relation [XT]. Whether w = —1 



or not will clearly signify whether the dark energy is constant in time or not. 
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FIG. 2. The reconstructed equation of state assuming Qm = 0.30 (solid curves). The curves 
assuming fijvf = 0.25, 0.26, . . . , 0.35 (from top to bottom) are shown by dotted ones. 

In Fig. 2, we show the reconstructed w(y) for cosmo logical models with constant equation 
of state w = 0,-1/3,-2/3,-1 to examine the effect of the ambiguity in fl M . The template 
fj is constructed by assuming Vt M = 0.30 for 0.60 < y < 1 (1.78 > z > 0). The template 
is fitted using the ansatz Eq.(|||) by minimizing J2 yi {^i ~ ^{Ui)) 2 'A? with N=10 data points. 
We also plot the reconstructed w(y) for Qm = 0.25,0.26, . . . ,0.35(from top to bottom) as 
dotted curves. We find that the uncertainty in Qm would enlarge an error in w(y) for small 
y (large z). In particular, when we overestimate VLm-, w(y) may diverge at some y where the 
denominator in Eq.(||) may vanish. However, an error in w remains relatively small near 
y — 1: 10% uncertainty in VLm results in at most 19% error in w for y > 0.80 [z < 0.56). The 
error is largest for the cosmological constant (19% error) and becomes less significant for 
larger w. For example, for w = —2/3 model, the error is less than 12% for y > 0.80. Hence 
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it might be possible to discriminate between w = —1 and w ^ — 1. The combination of SNIa 
measurements and high precision measurements of the power spectrum expected from the 
Microwave Anisotropy Probe (MAP) and Planck satellites could make a clear distinction 

El- 



B. reconstructing the effective potential 

Reconstructing equations are 

3 , , t , ^d 2 r(z)/dz 2 3 n n , , 




+ ( X + z ) tJ rT73 n3 - o^m(1 + *) 3 , (6) 



(dr(z)/dzY (dr(z)/dz) 6 2 



(dr(z)/dz) 

(i + zY 



. d 2 f(z)/dz 2 n „ 



(7) 



where = k0, V = k 2 V/H 2 with k 2 = 87rG. Alternatively, in terms of y = l/yl + 



z 



V( ) = 6 _ ^r/dy 2 _ 3Qm 
lyj y 6 (dr/dy) 2 y 5 (dr/dy) 3 2y 6 ' U 

M^/)\ 2 _ 12 | 4rf 2 f/A/ 2 3 ^ (rff/^) 2 
\ dy J y 2 y(dr/dy) M y 2 

In order to demonstrate the effectiveness of the fitting function Eq.(|3[), we reconstruct 
the effective potential of the quintessence field. We consider three kinds of potentials which 
have some theoretical backgrounds; (a) cosine type jllj, V(cj)) = M 4 (cos(0) + 1); (b) inverse 
power law type |TJ, V($) = M 4 0" a ; (c) exponential type |8|, V(</>) = M 4 exp(-A0). M 
and A are fixed to give Qx = 1 — Qm- ct = 4 is assumed hereafter. 

The results are shown in Fig. 3. The dotted curves are numerically reconstructed w(y) 
and V((p) with VLm = 0.30 being assumed, while the solid curves are the original ones up to 
y = 0.6. We fix the present value of (f> to unity. 

We also plot the dashed curves assuming Q M = 0.25,0.27,0.33,0.35. Since the smaller 
y, the larger the error in w(y), the range of significantly depends on the assumed VL M . We 
note that for Vt M > 0.30, w becomes less than —1 and thus the right hand side of Eq.(||) 
turns negative at some y. So we stop the reconstruction of V(<f) there. The range of <fi is 
larger (smaller) for smaller (larger) Qm- For example, in the case of cosine potential, the 
true range of is 0.41 < < 1.00, while the reconstructed range is 0.41 < < 1.00 for 
Q M = 0.30, 0.17 < < 1.00 for Q M = 0.25, 0.25 < < 1.00 for Q M = 0.27, 0.59 < } < 1.00 
for Q M = 0.33, and 0.66 < <p < 1.00 for Q M = 0.35. However, it is interesting that the whole 
shape of V(<f>) is less sensitive to the uncertainty in Qm, although w and the range of <f> are 
dependent on the assumed value of Qm- If we assume smaller Qm, the resulting potential 
energy is larger via Eq.([|). On the other hand, decreases more rapidly back in time via 
Eq.(||). The opposite is the case for larger Qm- Both effects make the reconstructed shape 
of V(4>) converge to the true one. 



6 




n — i — I — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



-0.3 




0.2 0.4 0.6 0.8 


Fig. 3(b) 




FIG. 3. The reconstructed equation of state and the effective potential: (a) cosine potential, 
(b) inverse power law potential, (c) exponential potential. Solid curves are the original equa- 
tions of state or effective potentials. Dotted curves are for VLm = 0.30. The curves assuming 
VLm = 0.25, 0.27, 0.33, 0.35 (from top to bottom for w{y) or from left to right for V) are shown by 
dashed ones. 



IV. RECONSTRUCTING THE EQUATION OF STATE AND THE EFFECTIVE 

POTENTIAL FROM SIMULATED DATA 

We simulate the actual reconstruction procedure using numerically generated data r\ = 
r {Vi) + $ r i (i — 1) • • • ) N) with 5ri being Gaussian distributed (zero mean and variance 
ar{i)i)). The simulated data assumes a cosmological model with Q M = 0.30. We consider 
N = 30 data and take a = 0.03 or a = 0.005. The former error is the distance error of the 
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binned data expected from observations of 200 supernovae by SNAP, while the latter is for 
6000 supernovae. We only consider statistical uncertainties. 

We distribute the data uniformly in y from y = 0.95 (z = 0.11) to y — 0.60 (z = 1.78). 
We perform thousands of Monte-Carlo realizations. The 68% confidence intervals of the 
reconstructed potentials are shown in Fig. 4 and Fig. 5. The horizontal axis is the averaged 
value of 4>, while the vertical axis is the averaged value of V. The effect of the reduced 
error may be dramatic. We allow for 10% error in Q]\j to assess the ambiguity in the 
reconstruction of the potential. Obervations of very distant supernovae at z > 3 by NGST 
(Next Generation Space Telescope) and/or observations of galaxy cluster abundance by SZE 
(Sunyaev-Zel'dovich Effect) survey will determine Om to a few % p0| . 

It should be noted that the fitted f (y) does not necessarily satisfy the positivity of the 
right-hand-side of Eq.(|5|), that is, the weak energy condition. We perform the reconstruction 
using only the data which satisfy the weak energy condition.^ That induces the bias toward 
larger 2 . This is why the range of is larger than that of the true one and why the inter- 
vals of the reconstructed potential is distributed downward. Such an effect is particularly 
significant for the cosine type potential because in this model the equation of state is almost 
w = — 1 at higher redshift. Therefore, good estimate of the upper bound of Qm is crucial 
for the success of the reconstruction. 



5 Conversely, Eq.(^) might provide an upper bound on Um, if we assume that the dark energy 
respects the weak energy condition, although it is not always the case |i~9| . 



S 



\ 1 



0.2 0.4 




-0.5 0.5 1 

Fig. 4(c) 

FIG. 4. One sigma intervals for the reconstructed quintessential potential assuming luminosity 
distance error of 3% with N = 30 data: (a) cosine potential, (b) inverse power law potential, (c) 
exponential potential. Solid curves are the original potentials. Dotted curves are for Q<m = 0.30, 
short dashed curves for Q,m = 0.27, long dashed curves for VLm = 0.33 
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V. SUMMARY 

We have studied the feasibility of reconstructing the equation of state of dark energy 
and the effective potential of the quintessence field from SNIa data by taking into account 
the uncertainty in Om as well as the error in the luminosity distance. We have found that 
w and the range of are dependent on the assumed value of Qm, while the whole shape of 
V{4>) is less sensitive to the uncertainty in Qm- If &>m could be constrained to some 10% 
accuracy by other observations, which may not be unrealistic expectation pD| , then future 
high precision measurements of distances to thousands of SNIa could reveal the shape of the 
quintessential potential. 
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